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We discuss the effect of tlie locai field polarization in metallic nanostructures where large field 
enhancements can be induced by excitation of localized surface plasmons. To demonstrate the 
importance of these effects, we study a few model systems in 2D and 3D. As an experimental probe 
of this effect, we suggest to use depolarization ratios in Surface Enhanced Raman Spectroscopy 
(SERS) measurements. We attempt to relate the polarization of the local field to the depolarization 
ratios in SERS conditions. This approach could also lead to some better understanding of the SERS 
mechanisms at a microscopic level. 

PACS numbers: 78.67.-n, 78.20.Bh, 78.67.Bf, 73.20.Mf 
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I. INTRODUCTION 

Raman depolarization ratios of molecular vibrations 
are one of the most fundamental properties of inelastic 
light scattering in liquids they are widely used in con- 
junction with other types of spectroscopies to infer indi- 
rectly the symmetry of molecular vibrational modes. Es- 
sentially, every vibration in a molecule can be assigned 
to a given irreducible representation (or linear combina- 
tions of them) of either the overall symmetry point group 
of the molecule or one of its sub-structures. Once the 
symmetry is known the Raman tensor of the vibration 
with respect to a fixed molecular coordinates system is 
also known. A collection of randomly oriented molecules 
will produce average intensities for parallel (/||) or per- 
pendicularly (/_l) polarized scattering configurations in 
Raman spectroscopy, from where a well defined depolar- 
ization ratio p = I±/I\\ follows. Both (/||) and (/_l) only 
depend on tensor invariants which do not change under 
rotations of the reference frame. This is a canonical text- 
book problem in non-resonant inelastic light scattering of 
liquids and poly-crystalline solids 1], which has also its 
exact counterpart in elastic (Rayleigh) scattering (with 
the linear polarizability tensor playing the role of the 
Raman tensor) j^] . There is ample theoretical and exper- 
imental understanding of this effect in Raman scattering, 
which started with the pioneering work of Porto, Skinner, 
and Nielsen^, Q in the 60's in the angular dependence 
of inelastic light scattering of simple liquids. 

Notwithstanding, under Surface Enhanced Raman 
Scattering (SERS) conditions in liquids (colloids) de- 
polarization ratios of Raman signals of the analyte 
molecules (typically dyes) are most of the time anoma- 
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lous and, in addition, show some laser wavelength depen- 
dence. R is interesting to understand the origin of these 
anomalies for such a fundamental physical property of 
the scattering process like the depolarization ratio. The 
problem is intimately related to the issue of local fields in 
SERS and in plasmon-supporting structures in general. 
SERS is made possible by highly localized plasmon res- 
onance excitations called hot-spots, which achieve large 
optical amplifications and boost Raman scattering sig- 
nals by several orders of magnitude^, |^. Highly (sub- 
wavelength) localized hot-spots have been proposed and 
studied both theoretically |3,I3 and experimentally 0,0, 
ITU . R is the purpose of this paper to shed some light on 
the polarization aspects of hot-spots and how they affect 
basic properties of the symmetry of the scattering. We 
also explore how these effects could be used in principle 
to learn more about the nature of hot-spots, local fields, 
and laser- field amplifications by plasmons. 

We would like to address in the paper two main issues: 

• Local-field polarization effects are often overlooked 
in most treatments of the SERS effect, but should 
be important for many essential properties of the 
scattering process. We shall show that depolariza- 
tion ratios are mainly determined by plasmons in 
a SERS environment. The most striking demon- 
stration of this effect is the breakdown of normal 
depolarization ratios in single isolated objects, as 
demonstrated in the next sections. 

• Depolarization ratios can, in principle, assert some- 
thing about the SERS enhancement mechanism it- 
self and, in addition, the equally important issue of 
molecular orientation and geometry of the analyte 
on the metallic surface. 

All these issues will be addressed with a variety of 
examples for different geometries and assumptions, to- 
gether with some experimental evidence for anomalous 
depolarization ratios at the end. 
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II. BASIC DEFINITIONS 

SERS enhancements by one or several objects depend 
on the geometry of the objects themselves, their symme- 
try, and their orientation in space with respect to the 
fixed scattering polarization directions in the laboratory 
frame. In addition, the scattering can come from a single 
analyte (molecule) in a specific place on the surface of 
the metal, or from the multiple contributions of a distri- 
bution of molecules in many places. Last, but not least, 
the Raman tensor of the analyte can be of many dif- 
ferent types (symmetries); hence the orientation of the 
main axes of the tensor with respect to the metal surface 
at specific points has to be specified also. In effect, the 
possibilities are endless. We need to start by imposing 
some restrictions on the type of objects, Raman tensors, 
and molecular orientations at the surface that we will 
examine. 

We first address the issue of the origin of the SERS sig- 
nal in terms of single/multiple dyes contributions. Inter- 
pretations of the data based on single molecule detection 
have been claimed many times in the past[3i [IB 13 ■ 
If single molecule detection is achievable then the po- 
larization properties of the local field at the position of 
the molecule will be the only parameter that counts. 
The same situation happens if the presence of analyte 
is widespread over the metal but there are high inten- 
sity hot-spots dominating the signal. In such a situation, 
the polarization properties of a specific place in space 
and/or the orientation of the Raman tensor of one or a 
few molecules are the main factors determining the sym- 
metry of the scattering event. The polarization proper- 
ties would then be an average (over time or space) of 
many of these hot-spots. Polarization studies in SERS in 
single molecule limits are very rare^^0| and not fully 
understood in our opinion. 

The most common experimental situation, however, is 
one in which there are many contributions from many an- 
alytes which are spatially distributed in different places 
and, accordingly, are exposed to different local fields. 
This is for example the case of colloidal solutions of gold 
or silver particles at relatively high (~ 100 nM) concen- 
tration of analyte and measured in immersion. If con- 
ditions are chosen so that the solution is stable, we can 
rule out the presence of large clusters of particles (which 
would precipitate) and the signal must therefore come 
from single particles or possibly pairs of particles (either 
stable pairs or dynamical pairs formed during close en- 
counters of two particles). In order to observe a good 
SERS signal from such solutions, higher analyte concen- 
trations are required, corresponding typically to ~10^- 
IC* molecules per particle. To study the polarization 
effects in such systems, one has to consider the signals 
of all these molecules and also average over all possible 
orientations of the particles. This is only possible in sim- 
ple cases and it is the main reason for the restrictions we 
shall impose in the modelling. 



A. Basic geometries and Raman tensors 

We shall discuss here a few model cases: spherical par- 
ticles, elliptical particles, and dimers (two particles in 
close proximity). In each case, we will model the SERS 
depolarization ratio, which can be measured experimen- 
tally, under different assumptions. 

There are three main factors that can affect the calcu- 
lated depolarization ratio p. 

• The Raman tensor of the probe molecules: We will 
consider the cases of an isotropic Raman tensor 
(/so) and of a uniaxial tensor. These are two ex- 
treme cases of a variety of molecular symmetries 
present in real analytes. Many common SERS an- 
alytes (like rhodamine 6G) have strongly uniaxial 
Raman tensors for most of their vibrations. 

• The adsorption configuration of the molecule: This 
is obviously an important factor when the Raman 
tensor is not isotropic. For the uniaxial case, we 
therefore need a further assumption about the way 
the molecule adsorbs on the metal surface and, in 
particular, how the main axis is oriented with re- 
spect to the surface. We will consider the two ex- 
treme cases of a molecule adsorbed with its axis 
tangential (Tan) or normal (Norm) to the surface. 
We will also consider the case where the adsorption 
geometry is random (Rand). In this latter case, all 
possible orientations of the molecules are consid- 
ered and averaged. 

• Finally, the way the molecule couples to the local 
field polarization: This is also a very important 
aspect because it is directly related to the SERS 
enhancements mechanism and has in our opinion 
been overlooked. 

We shall have a small digression on this latter issue 
here. If we consider a molecule with a highly uniaxial 
Raman tensor oriented along e„i, then the excited dipole 
has the form (case A): 

jf- ^ a (^E ■ er,^^ Crn, (1) 

where E is the local electric field and a is the main com- 
ponent of the uniaxial polarizability. This is a case of- 
ten found in normal Raman scattering. The intensity 
scattered and detected in the far field in a given direc- 
tion Cd, along a given polarization Cp (cp _L e^) is then 
/ cx |p • Cpl^. However, if we use this expression in SERS 
conditions, the SERS signal is only proportional to the 
power of two of the local field enhancement factor (and 
not power of four as usually assumed(s|). In order to ob- 
tain a power of four in the scattering efficiency, we need 
to model the enhancement process due to re-emission of 
this dipole. One could use the following expression (case 
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B): 



and 



a\E\ 



E ■ 



(2) 



Anywise, this implies that the enhancement in re- 
emission has no polarization dependence. Another view 
is to assume that the excited dipole couples only to the 
local field enhancement along its direction; we then get 
(case C): 



(3) 



In both cases, we have also assumed that the local field 
enhancements relevant to re-emission were the same as 
those derived for excitation. This is a common assump- 
tion but is not necessarily valid. Assuming we know the 
local field enhancements for a given geometry, we can 
then calculate the polarization property of the SERS sig- 
nal for each of the cases A, B, and C at each point of the 
surface S of our metallic structure/s. These can then be 
integrated over the surface, assuming a continuous dis- 
tribution of probe molecules of type Iso, Tan, Norm, or 
Rand. 



B. Symmetries and averages 

As an example, and to provide a basis for further dis- 
cussion, we carried out calculations for simple model sys- 
tems in 2 dimensions (2D). The structures are in the 
plane (yOz) and the exciting beam is coming from the 
perpendicular direction (Ox), with a field polarization 
along z. The local fields at the surface are calculated 
within the electrostatics approximation. The calcula- 
tions were carried out numerically using finite element 
modelling described elsewhere po|. We assume that the 
signal is detected in the far field along {Ox) and analyzed 
for parallel {Oz) and perpendicular {Oy) polarizations. 
For a given structure, we calculated /|| and for each 
of our assumptions, ft (i) denotes the unit normal (tan- 
gential) vector at any point on the metallic surface. We 
therefore used the following expressions: 



Iso— A 



(4) 



jN orm — C 

'ii 



E ■ n 



nIdS. 



(8) 



Similar expressions for Tan— A, —B, —C replacing n with 
t are also obtained. Expressions for I± are obtained by 
replacing z subscripts with y. We will also for each case 
estimate the total SERS intensity by means of: 



Irot — I\ 



T 



(9) 



For structures without rotational symmetry, the cal- 
culations were repeated for 50 different orientations, pa- 
rameterized by an angle (3. The total signals are then 
obtained by summing over all possible orientations. In 
our structures, we need only to consider f3 from to 7r/2 
by symmetry, so we have: 
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7r/2 



h{P)df3, 



(10) 



and the depolarization ratios is then calculated according 
to 



(11) 



In most simulations of metallic structures presented in 
the literature, such as for ellipses or dimers, only one 
or two possible orientations are considered. In liquids, 
contributions from all orientations are averaged and the 
resulting intensity profile can be very different to that 
of a single configuration. Independent of polarization 
effects, we believe that the orientation-averaged SERS 
intensity profiles presented here provides a new insight 
into plasmon resonances in liquids. 

For the Rand cases, one has to consider all possible 
molecular orientations. In 2D we restrict ourselves to: 



cos(7)e*y -I- sin(7)4 



(12) 



Estimating the intensities and averaging over 7, one finds 
that the results Rand — B and Rand — C can be directly 
obtained from that of Iso ~ B (also equal to Iso — C) , 
namely: 



^Rand—B 



Iso-B 



1 



„Iso-B 



(13) 



jIso—B jIso—C 





E 







\E,rdS, 



(5) 



jRand—B 
^Tot 



- jIso-B 
i^Tot ' 



(14) 



jN orm — A 

^11 



E ■ n 



n'^dS, 



(6) 



and 



Rand—C 



bp 



Iso-B 



„Iso-l 



(15) 



rNorm — B 







1 




L 


E 




E-n 











n^dS, 



(7) 



jRand-C ^ rIso-B 



'-Tot 



'^Tot 



(16) 
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As a reminder, in non-SERS conditions in 2 dimen- 
sions, p should be for an isotropic molecule and 1/3 for 
the uniaxial case {Rand). 

For all calculations, we used the bulk dielectric con- 
stant of Ag (a typical metal for SERS applications) as 
reproduced by a Drude model with the best fit to tabu- 
lated data in the region of interest [23, namely: 



Jso-B 



|fc|V2+|fc|V2 



1 + |/s|4/2 + 3|fc|2/2 + 2Re(fc)2 



Norm — A Tan— A 



Norm-C Tan-C 



1/3, 



1/5, 



(22) 



(23) 



(24) 



e(A) 



1 



(17) 



N orm— B 



where too = 4, Ap = 141 nm, F = 17000 nm, and A is 
the wavelength. The surrounding medium is water with 
em = 1.77. 



III. THE SIMPLEST CASE: SERS BY A SINGLE 
SYMMETRIC OBJECT 

A simple object like a disk in 2D or a sphere in 3D 
has a well defined main plasmon resonance. If the size 
of the object (a) is much smaller than the wavelength 
(a <C A) of the light, then the electrostatic approxima- 
tion is valid and the resonance occurs at the frequency 
where Re(e(A)) « — e™, or — 2em for a disk in 2D or a 
sphere in 3D, respectively. If (a ~ A) on the contrary, re- 
tardation effects and higher order multipoles play an in- 
creasing role and define additional resonances. The cur- 
rent understanding of the SERS effect is that it benefits 
from red-shifted plasmon resonances which are either 
shape/size- related or, more often than not, coming from 
plasmon-plasmon interactions among particles|l3l Il4| . 
These latter resonances are believed to be the dominant 
contribution to the SERS signal in general. 

The case of a 2D sphere (or disk) of radius a, is partic- 
ularly easy because of the rotational symmetry (no need 
for different /3's in the averaging). It can also be eas- 
ily solved analytically in the electrostatic approximation. 
We consider an Ag-disk with complex dielectric function 
e(A) in water (e^ = 1.77). The important parameter in 
such problem is the relative dielectric constant tr = e/e„i. 
For an incident field Eq = le^, the local field at a point 
on the surface with coordinates a, 9 is simply the super- 
position of the incident field and of an induced dipolar 
field: 



EyiO) = fcsin(26l), 
E^{e) = 1 - fccos(26'), 



with 



k 



Cr — 1 
Cr + 1' 



(18) 
(19) 



(20) 



We can therefore derive the following analytical expres- 
sions: 

Iso-A_ 1^172 



and 



P 



5|eJ2 + i 



5' 



(25) 



(26) 



P 



l + |A:|2/2' 



(21) 



where pR'^™^-B ^-^^^ pRand-c ^g^^^ g^jg^ -^^ evaluated an- 
alytically from the expression of p^'^"^^ using Eqs. H13|l 
and ifTCjl . 

FigureQlshows the wavelength dependence of the depo- 
larization ratio under the various assumptions considered 
here. An important aspect to help distinguish these cases 
is also the total SERS intensity, which is not reflected in 
the values of p. We know that the SERS signal is usually 
strong enough to be detected, and any valid assumption 
should therefore predict a reasonably strong SERS signal 
in addition to a correct depolarization ratio. We show in 
Fig. n the calculated SERS intensity as a function of 
wavelength. As can be expected, the three cases where 
the SERS enhancement is only proportional to the square 
of the field amplitude (cases A) show a relatively small 
SERS intensity. These three cases will therefore not be 
studied any further. All the other cases show a strong en- 
hancement at the localized surface plasmon resonance of 
the object (Re(e) — — Cm, A ~ 339 nm in this case). For 
Tan — B and Tan — C , this enhancement drops sharply 
beyond this resonance, while a small enhancement still 
exists in the other cases. 

This is a very interesting case for more than one rea- 
son. It is interesting to see how a large variety of cases 
is observed for p in an object which is arguably one of 
the simplest SERS enhancing object we can think of. In 
some cases, p — \/Z and is independent of wavelength 
{Norm - A, Tan - A) or 1/5 {Norm - C, Tan -C). In 
other cases, it increases to 1 close to the surface plasmon 
resonance and decreases slowly beyond that, towards a 
limiting value of 1/5 {Iso-B), 5/13 {Rand-C), or 1/2 
{Rand — B). Tan — B and Norm — B show a somewhat 
different behavior, the first showing an increase towards 
1, while the latter peaks sharply at 1 before the surface 
plasmon resonance. 

Conceptually, it is interesting to see that for the case 
of an isotropic probe (which should exhibit p = in 
free space) p can be as large as 1 as a result of the 
interaction with the metallic disk {Iso — B aX reso- 
nance). This is surprising because the disk itself is also 
isotropic, but the anisotropy is introduced by excitation 
of non-isotropic dipolar plasmon resonances in the object. 
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This demonstrates clearly that the localized plasmon res- 
onances are crucial in the polarization mechanisms and 
cannot be excluded from the discussion. It is even clear 
from Fig. ^ that the plasmon resonances, and not the 
probe molecules, are the dominant factor in the origin of 
the depolarization ratios; what counts is the symmetry of 
the plasmon resonance for the polarization properties of 
the scattering. This summarizes in many ways one of the 
most important points of this paper. 

It is interesting to see, in addition, the reasons for the 
influence of plasmon resonances in terms of field patterns 
for the simplest possible case treated in this section. In a 
way, all the other effects are further degrees of sophisti- 
cation of this simple example. Figure |21 shows for this 2D 
case the intensity patterns on the surface for an incident 

polarization along z of \Ez\^, {Eyf^ , and the total field 
^ 2 

E , at three different wavelengths chosen to be below, 
at, and above the main surface plasmon resonance of the 
object (A = 339 nm). The case can be solved analytically 
and all the patterns can be interpreted as a direct super- 
position of the external field and an equivalent induced 
dipole on the disk. It is the complex polarizability of the 
metal that makes the dipole response non-trivial because 
of the relative amplitude and phase shift of the response 
with respect to the exciting field. From Fig. |21we can see 
that even if we had molecules with an isotropic Raman 
tensor (emitting in the same direction as the excitation), 
the depolarization ratio will be changing with wavelength. 
Ignoring the Stokes shift, the integrated intensity on the 
surface for \Ey\^ divided by that for \Ez\^ will give the 
depolarization ratio at that wavelength. It is easy to infer 
by visual inspection that the depolarization ratio will be 
small at A = 312 nm, large at resonance (A — 339 nm), 
and small again at longer wavelengths (A = 364 nm). 
It is large (close to 1) at resonance because the induced 
dipolar response completely dominates the signal. The 
(analytical) example in Fig. [3 is possibly the simplest 
example of plasmons affecting the Raman depolarization 
ratio by changing the nature of the coupling between the 
laser and the molecules through the local field. In prac- 
tise the Raman tensor of the molecule itself will play a 
role, adding an additional degree of complexity to the 
problem. 

As we shall show later for a particular analyte with 
a uniaxial Raman tensor, SERS experiments on col- 
loidal solutions show that the depolarization ratio indeed 
changes with wavelength and ranges from 0.4 to 0.7 in our 
measurements. This observation therefore rules out some 
of the cases studied here. We can identify three cases 
which are compatible with this observation and predict 
a reasonable SERS intensity: Iso — B, Rand — B, and 
Rand — C. In order to determine which of these cases is 
the most realistic, one first needs to look at more real- 
istic models. It is in particular important to determine 
the effect on p of different particle shapes and to assess 
the validity of using a 2D model to infer 3D properties. 
Finally retardation effects (not taken into account in the 



electrostatics approximation) can have some influence on 
the plasmon resonances and therefore also on p. We ad- 
dress some of these issues in the next sections. 



IV. SHAPE EFFECTS 

To understand at least qualitatively the effects of 
shape, we will consider the case of an ellipse with an 
aspect ratio between long and short axes of 1.5. Be- 
cause we now loose the rotational symmetry, we consider 
the contributions to the SERS signal from ellipses with 
50 possible orientations. Figure 13 shows the calculated 
total SERS intensity under different assumptions. As ex- 
pected, the plasmon resonance is now split into two reso- 
nances: one red-shifted at 363 nm, the other blue-shifted 
at 320 nm. In many instances in the literature, excita- 
tions along one or the other principal axis are usually 
considered separately and only one resonance is observed 
at a time. Both resonances are observed here because all 
possible orientations are considered and summed. Figure 
13 also shows the predicted depolarization ratios. Again, 
a wide variety of different behaviors are observed, but we 
will concentrate on the three cases identified previously: 
Iso — B, Raiid — B, and Rand — C. All three show a 
similar pattern with p peaking for A = 340 nm to a value 
at 0.8 - 0.9, then at A = 388 nm to a value of 0.6 - 0.8, 
and then decreasing to the same limiting value as that 
observed for the disk. 

It is interesting to note that the peaks of p do not cor- 
relate exactly with the plasmon resonances observed in 
the SERS intensity. The first peak remains at the plas- 
mon resonance of a sphere, while the second appears at 
even longer wavelength than the red-shifted plasmon res- 
onance. Also, the maximum value of p is now slightly 
smaller than 1. This can have some important conse- 
quences for colloidal solutions with a distribution of sizes 
and shapes. For a given excitation wavelength Al, the 
SERS signal should have a larger contribution from parti- 
cles in resonance with this wavelength. If there are many 
such particles, the depolarization ratio will be determined 
by these and will then correspond to the dip in Fig. 13 at 
363 nm (i.e. p is smaller than its maximum value). How- 
ever, if Ai is longer than most particle resonances, the 
measured p should correspond to values to the right of 
the dip in Fig. |3| This could lead to an increase in p if 
Al remains close to most resonances, or to a decrease at 
much longer wavelength. Such considerations are impor- 
tant when interpreting experimental results in colloidal 
solutions because polydiversity is inevitable. 



V. CASE OF A 3D SPHERE 

It could be argued that conclusions gained from 2D 
models on depolarization ratios should be taken with 
care; it is necessary to estimate the changes that will 
arise from looking at real 3D cases to gain some insight 
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on how reliable the qualitative predictions of 2D models 
are. Numerical simulations in 3D are much more com- 
putationally demanding than the 2D cases. To assess 
the importance of the corrections arising from a full 3D 
model, we will therefore study the case of a 3D sphere and 
compare qualitatively the results to the 2D case. Simi- 
lar to the 2D disk, the 3D sphere problem can be solved 
analytically in the electrostatics approximation. For an 
incident field Eq = le^, the local field outside the sphere 
is simply the superposition of the incident field and of an 
induced dipolar field: 



E = 3KCOs{9)er + (1 — njcz, 



(27) 



where (r, 9, (jy) are the standard spherical coordinates and 



(28) 



Calculations can be carried out in a similar way as for 
the 2D case. To illustrate this, we describe below the 
Norm — C, Tan — C, and Iso — B cases. 

In the Norm — C case, the molecules are assumed to 
be adsorbed with their main axis normal to the metal 
surface. The parallel and perpendicular intensities are 
therefore given by: 



/,f = I \E- n\*{n ■ e.fdS = ^ |1 + 2Kp , (29) 



and 



lNor,n-c ^ / \E.n\'^{n.eyfdS = — \l + 2K\''. (30) 
Js 35 

The depolarization ratio is then independent of wave- 
length and equal to pN"™-c — 1/5^ which is the same 
result as that obtained in 2D. 

For the case of Tan — C, we assume the molecule is 
adsorbed with its axis tangential to the metallic surface. 
However, in 3D, this does not define completely the ori- 
entation of the axis, and we therefore need to average 
over all possible orientations of the molecular axis (tan- 
gential to the surface). For this, we assume the axis to 
be along the unit vector: 



e,n = cos7ee -|- sin 76*0. 
We then calculate: 



(31) 



jTan-C 
1 



(7) = / \E-e,n\\em-S,YdS (32) 
Js 

= -cosSll-'^l', (33) 



and 



(7) = / \E-e„,\\e„,-eyfdS (34) 



|l-K|''cos*7(l-H6sin^7). (35) 



105 



Averaging over 7 (0 < 7 < 27r), we obtain: 



jTan-C ^ 



and 



jTan-C 



35 



(36) 



(37) 



The depolarization ratio is therefore again p2"an-c _ 
1/5, independent of wavelength and equal to the 2D case. 

In the 2D section, we identified the cases Iso — B, 
Rarid — B, and Rand — C as the most interesting. The 
Iso — B case can be calculated analytically, to wit: 



Iso-B 



II 105 

[(3|e^|2 -t- 4)(5|e^|2 -2 + 4Re(eO) + 56] 



(38) 



and 



I'r''' = - ^ (39) 



Hence, the ratio: 



Jso-B 



P 



|e,-l|2(3|e,|2+4) 



(3|e^|2 + 4)(5|e^|2 - 2 + 4Re(e,.)) + 56 ' 



(40) 



For the Rand — B and Rand — C cases, one needs 
to consider molecules with a random orientation of their 
axis: 

e„i = sin7Cos5e2: -I- sin 7 sin (5 ej, -1-003762, (41) 

with < 7 < TT and < (5 < 2ti. Averaging over these 
parameters, one can show that: 



Rand— B 



Iso-B 



3 -t- 2pi^°-B ' 



rRand—B 



12 
15' 



Iso-l 



^Rand—C 



l + 6p 



15^11 



Iso-B 



Iso-B 



2p 



Jso-B ' 



and 



~ 35^^ 



15 



Iso-B 



(42) 



(43) 



(44) 



(45) 



These results are plotted in Fig. ^ with the equiva- 
lent 2D results for comparison. We first notice, as ex- 
pected, that the resonance profile of the SERS intensity 
is further red-shifted in the 3D case, peaking at 388 nm 
(Re(e^) = -2) instead of 339 nm (Re(er) = -1). Re- 
garding the depolarization ratios, they show the same 
qualitative behavior as in 2D with some small quanti- 
tative changes: they are smaller for the 3D case below 
the plasmon resonance, and slightly larger beyond this 
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resonance tending towards the same value at long wave- 
length. Also, the maximum depolarization ratio is no 
longer 1 as in the 2D case but around 0.7-0.9 depend- 
ing on the case. Finally, the maximum of p now occurs 
slightly before the resonance at around 368 nm. For the 
cases Rand — B and Rand — C , p remains quite high, at 
least larger than 0.5 beyond the resonance, while it goes 
down to around 0.3 for Iso — B. 

The overall conclusion of translating the simulation 
from 2D to 3D shapes is that the essential qualitative 
features are already contained in the 2D simulation. The 
quantitative picture will, of course, be different. But the 
core qualitative (and semi-quantitative) effect of the main 
plasmon resonance is already contained in the 2D simu- 
lation. 



VI. RETARDATION EFFECTS 

So far, all the calculations have been performed within 
the electrostatics approximation. This is strictly speak- 
ing valid for very small objects a <C A, where the electric 
field can be assumed to be constant over the object. In 
practise, retardation effects can become non-negligible 
in the visible range for objects of dimensions « 30 nm 
(A/ 10). Retardation effects tend to red-shift the plas- 
mon resonances and add more structure to the resonance 
profile. In order to appraise the effect of retardation on 
the depolarization ratios, we used Mie theory .22] to cal- 
culate exact results for the simple case of 3D spheres. 
We follow the treatment of Mie theory given by Bohren 
and Huffman [2^. For a sphere of radius a, we calcu- 
late the exact value of the field at the surface E{a, 9, (f). 
We can then evaluate the relevant integrals to calculate 
SERS intensities and depolarization ratios as before. In 
particular in the Iso — B case, we calculate: 



and 



tIso — I 



\E\''\EA-'dS, 



- / \E\^\Ey\''dS, 



Ji--^= / \E\''\EJ''dS. 



(46) 
(47) 
(48) 



/j_2 is a magnitude that cannot be measured when de- 
tecting along the excitation direction x, but will be useful 
nevertheless to calculate other ratios. The depolarization 
ratio is then 

jIso-B 

= yra- (49) 

By considering molecules of random orientations, one can 
also show that: 

2 



rRand— B 



rRand— B 



15 
2 

15 



Q rIso—B , jIso — B , jIso — B 
jlso—B I o jlso—B I jIso — B 



, (50) 

, (51) 
(52) 



fRand—C 



jRand—C 



35 
2 

35 



^Jlso-B ^ jIso-B ^ jIso-B 



tIso- 



r jIso — B , tIso — B 



(53) 

(54) 
(55) 



The corresponding depolarization ratios p^i^^^d-B ^^^^ 
pRand-c thcu be evaluated. 

In Fig. |31 we display the results for the Iso — B case, 
for sphere radii in the range a = 1 to a = 60 nm. As 
expected, the results for a = 1 nm are identical to those 
obtained from the electrostatics approximation. As a 
increases, retardation effects lead to a number of modifi- 
cations: First the plasmon resonance is split into several 
resonances, one of which red-shifts more the larger the ra- 
dius. The maximum SERS intensity also decreases as a 
increases. Concerning the depolarization ratios, the most 
interesting point is that it shows virtually no change be- 
yond its original resonance at 368 nm. In particular, the 
most red-shifted resonances of the intensity (which are 
the ones most relevant to SERS) do not affect p. Instead, 
p is modified in the regions of the secondary resonances. 
A closer look at the data indicates that p shows a dip for 
each of the resonances in the intensity (except the most 
red-shifted one). For example, for a = 30 nm, two dips at 
358 and 374 nm are observed for p, and they correspond 
to peaks in the intensity profile. However, the most red- 
shifted resonance (at ~432 nm) does not affect the profile 
of p. We conclude that in the regions of interest to SERS, 
retardation effects do not in general affect the main con- 
clusions on depolarization ratios gained from the simpler 
electrostatic approach. 



VII. COUPLED RESONANCES 
DIMER 



CASE OF A 



Finally, it is widely believed that the largest SERS en- 
hancements originate from interactions between two or 
more particles, leading to coupled plasmon resonances. 
In many situations, it is likely that the SERS signal origi- 
nates mainly from such interactions. We would therefore 
like to enquire whether such coupled resonances could 
affect our conclusions on the depolarization ratios. To 
do so, we calculated the depolarization ratio in the elec- 
trostatics approximation for two closely spaced 2D disks 
separated by d = 0.1a (where a is the radius). The re- 
sults are summarized in Fig. El The structure of the 
SERS intensity profile in the region of the single plasmon 
resonance (from 300 to 380 nm) is very complex. The 
coupled plasmon resonance is substantially red-shifted to 
440 nm in this case (404 nm for a separation oi d = 0.2a, 
not shown). This resonance is also the one exhibiting 
the highest integrated SERS intensity, as expected. The 
complexity of the intensity profile is due to the averaging 
of all possible orientations. In the literature, no orienta- 
tion averaging is usually considered, hence the intensity 
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profiles exhibit much less structure in general. A wide 
variety of depolarization ratios is observed depending of 
the assumption made. Some of the discarded models in 
this section (like Tan — A, Tan — B, and Tan — C) pre- 
dict values for p well above 1, a situation not observed 
experimentally. The most surprising result of the inclu- 
sion of coupled plasmons in dimers is that they seem to 
affect only the SERS intensity but not so much p, as can 
be observed in Fig. When compared to the results 
for a single disk (see Fig. [7| for example), coupled reso- 
nances appear to narrow the resonance profile of p and 
make it tend faster towards its limiting long wavelength 
value. This leads to a small decrease in p for intermediate 
wavelengths. 

Some of the other models not shown in the figure like 
Norm — B and Norm — C, show a p that tends quickly 
towards its limiting value of 1/5, which was also predicted 
for disks and spheres (2D and 3D) and ellipses (2D). A 
value oi p = 1/5 should therefore be observed in most 
SERS conditions if one of these assumptions were correct. 
This is not the case, for experiments with a particular 
analyte like BTZ2 (used in the experimental section later 
on), but could be correct in other cases. We therefore 
conclude that the assumptions Norm — B, and Norm — C 
for the interaction of the dye with the metal cannot apply 
for the specific analyte we shall explore later. 



VIII. SOME EXPERIMENTAL EVIDENCE 

We have already ruled out a number of possibilities 
and now concentrate on the most likely ones in terms of 
a possible connection with experimental results, namely: 
Iso — B, Rand — B and Rand~C models. All three cases 
show a qualitatively similar behavior, which is summa- 
rized in Fig. [3 for the case Rand — C for various model 
geometries. We can draw from the previous calculations 
a number of conclusions about the depolarization ratios 
in these cases: 

• In a SERS experiment on colloidal solutions 
with non-interacting colloids, where single colloids 
would contribute to a substantial fraction of the in- 
tensity, one should in principle observe signals from 
a distribution of spheroidal colloids, and possibly 
some elongated ones. One therefore expects to ob- 
serve a p of a maximum of about 0.7-0.9. How- 
ever, we saw in two instances (2D ellipses, and 
3D spheres with retardation), that the maxima 
of p corresponded to dips in the SERS intensity, 
and vice-versa. If the excitation wavelength Xl is 
within the range of particle resonances, resonant 
particles will contribute more to the signal, and a 
smaller value for p should be observed (we are in 
the dip associated with a maximum of SERS inten- 
sity). When Al is increased beyond most particles 
resonances, p should go through a maximum and 
eventually decrease slowly towards a limiting value 



depending on the model chosen. However, this lim- 
iting value is very difficult to measure, because the 
SERS intensity also decreases substantially. 

• In interacting colloid solutions, SERS signals com- 
ing from interacting colloids will overwhelm that 
coming from single particles. Since the predicted 
depolarization ratio of a dimer vary little with 
wavelength at the coupled red-shifted resonance 
and beyond, one would expect a constant p at 
long wavelengths, p should also increase as Al is 
decreased towards the single (uncoupled) plasmon 
resonance. Other scenarios could be considered: 
for example, because the signal comes from a col- 
lection of various interacting objects, some interac- 
tions may be more dominant at some wavelengths. 
In a simple dimer picture, one expects that as the 
wavelength increases, more coupled (and therefore 
more red-shifted) resonances are dominating the 
signal. A calculation of 2D dimers with different 
separations shows that increased coupling should 
have a negligible effect on p. Another possibility is 
that the colloids may be able to merge and form 
a single elongated particle. Their behavior could 
then be more similar to that of an ellipsoidal par- 
ticle. 

Far from attempting a full experimental study of the 
effects reported here for reasons of space, we give an ex- 
perimental indication of some of the important aspects 
that can be observed experimentally which do support 
the idea that depolarization ratios in SERS are affected 
by plasmon resonances. Accordingly, we measured de- 
polarization ratios of SERS active liquids under various 
conditions and at different wavelengths. 

The experimental details are identical to those re- 
ported in Ref. [2^ , except for the analyte which in these 
experiments is the newly developed benzotriazole dye: 4- 
(5'-azobenzotriazolyl)naphthalen-l-ylamine (BTZ2)'265. 
The samples were prepared by mixing 0.5 niL of silver 
colloids 25] with 0.5 mL of water or KCl (concentration 
10 mM), and adding 40 of 5 BTZ2 to obtain an 
end concentration of 0.2 /iM of dye. The samples were 
then diluted by a factor of 4 (in water or KCl) to test for 
a possible effect of with multiple scattering on the depo- 
larization ratio. No evidence for this were found when 
comparing diluted and non-diluted samples. To calcu- 
late depolarization ratios, the signals were averaged over 
periods between 5 and 20 minutes, to avoid problems 
associated with SERS fluctuations. The ratios is then 
estimated for 3 consecutive experiments to check for re- 
producibility and estimate errors. 

These samples are stable for several weeks at room 
temperature. There is substantial experimental evidence 
of poly-dispersion in such Ag colloidal solution. Figure |S| 
gives a more definite idea of the complexity of real col- 
loids, like the ones we used in the experiments reported 
below, in terms of size and/or shape distribution. The 
colloids have an average radius of a = 30 nm and exhibit 
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an absorption spectrum peaking at 454 nm. Spheres of 
a — 30 nm should have a resonance at 430 nm, and the 
observed red-shift in absorption is therefore attributed 
to variations in the shape of the particles. Without KCl, 
there is a strong coulombic repulsion among particles, 
and the SERS signal originates from single particle reso- 
nances. The addition of KCl results in a screening of this 
repulsive force, making the occurrence of closely spaced 
pairs more likely either by collisions or direct formation 
of dimers or small clusters. This is enough to observe a 
large increase in SERS signal, which shows that interact- 
ing particles dominate the signal. Such samples remain 
however stable with no aggregation into large clusters. 



BTZ2 has several Raman peaks, the most prominent 
being around 1410 cm~^. We measured the depolariza- 
tion of this peak for a high dye concentration in water (no 
colloid) with 514 nm excitation and obtained p = 0.33, 
indicating a uniaxial Raman tensor. The depolarization 
ratios for this peak measured under SERS conditions are 
summarized in Fig. |3 as a function of laser wavelength, 
for samples without KCl (single particles) and with KCl 
(interacting particles). It is interesting to note that un- 
der SERS conditions, the depolarization ratios of all the 
other peaks was nearly the same as this one; a situation 
not happening in the pure dye and confirming that p is 
mainly determined by plasmon resonances rather than 
by the Raman tensor of the peaks. The background (or 
continuum, a well-known characteristic of SERS) also ex- 
hibits the same depolarization ratio, indicating its ori- 
gin is strongly related to that of the SERS signal. Our 
experimental findings are consistent with the previous 
discussion for this specific experimental implementation. 
The models Rand — C, Rand — B, or Iso — B show a 
good qualitative agreement. Rand — C shows the best 
quantitative agreement, while Rand — B seem to be pre- 
dict too high a depolarization ratio at longer wavelength. 
Iso — B is not relevant to this specific experiments, since 
we know from the depolarization ratio measurement in 
the pure dye that the tensor is not isotropic, at least when 
not adsorbed on the metal. Note the clear presence in 
Fig. I^of a dependence on both the laser excitation and 
the KCl concentration, which reveals the role of plasmon 
resonances brought about by interactions among colloids. 
Both effects are expected from the general considerations 
presented in this paper. We discuss briefly the results in 
the next section. 



We believe that further experiments on depolarization 
ratios could prove useful to the general understanding 
of SERS. In particular, experiments where parallel and 
perpendicular polarization intensities can be monitored 
simultaneously would enable one to study the fluctua- 
tions of p with time. This could provide insight into the 
origin of the SERS fluctuations, and a test of their inter- 
pretation as coming from single molecules. 



IX. DISCUSSION AND CONCLUSIONS 

The rotationally invariant components of the Raman 
tensor measured by depolarization ratios have a limited 
amount of information about the symmetry of the vibra- 
tion. If a new molecule is measured in liquid, of which 
no previous information is known, we can infer the sym- 
metry of the Raman tensor only to a certain extent. In 
general a combination of techniques are required. The 
same limitations occur for the changes in depolarization 
ratios by the presence of plasmon resonances: there is 
only a limited amount of information we can infer about 
the orientation of the molecule on the surface and it is 
very likely that we cannot decide between two competing 
models in certain cases. Like in the original case without 
plasmons, a combination of techniques will be needed for 
that task. 

Still, what this paper has unequivocally shown is that 
plasmon resonances produce "anomalous" depolarization 
ratios in standard SERS conditions and that the presence 
of the plasmons cannot be neglected or ignored. The key 
issue highlighted in this paper can be summarized with 
the textbook example of local fields presented in Fig. |21 
SERS is about local fields, and the polarization of the lo- 
cal field does not follow directly from that of the incident 
field. In that respect, we can safely say that plasmons 
are the most important factors controlling the depolariza- 
tion of the scattering process under SERS conditions and 
are responsible for the large number of anomalies seen in 
experiments; a large fraction of which have not been ex- 
plained or remain unaccounted for. The experimental 
results of the previous section do suggest that both the 
excitation wavelength and the state of aggregation of the 
colloids (producing coupled plasmon resonances) play a 
role in the depolarization ratios that are observed exper- 
imentally. 

In this respect, we also believe that the interpretation 
of depolarization ratios put forward in several previous 
publications 0,^3 will probably have to be revised. De- 
polarization ratios have been used in the past to infer 
something about the orientation of the Raman tensor of 
a few (or one) molecules, completely ignoring the fact 
that the laser and Stokes field couple to the local field 
through the local plasmons and their symmetries. We 
believe that a re-interpretation of the data including the 
role of the plasmons might be necessary in those cases. 

Finally, the limited amount of experimental informa- 
tion we have at the moment seems to to favor model C 
with respect to B for the scattering efficiency. This is a 
longstanding unresolved issue in SERS: whether the en- 
hancement operates separately for the laser and Stokes- 



field thus producing a factor of cx 



E 



in the scattering 



efficiency, or whether it is only the enhancement at the 

2 



laser frequency cx 



E 



followed by a re-emission process. 

Far from resolving this issue here, we suggest that de- 
polarization ratios could help to discard possibilities and 
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learn a bit more about the SERS mechanism at a micro- 
scopic level. This is a unique aspect of SERS provided 
by the fact that the probe is in direct contact with the 
local field. 
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FIG. 1: Depolarization ratios and total SERS intensities cal- 
culated for a 2D disk from different assumptions as a function 
of wavelength. Tan, Norm, Iso and Rand define the type of 
Raman tensor under consideration, as specified in the text 
(Sec. II), while A, B, and C account for the different models 
of the scattering process defined in Eqs. Q, (|5J, and 
Note that the SERS enhancements are plotted on a log-scale. 




FIG. 2: Intensity patterns in a 2D disk at three different 
wavelengths (columns) for the field component along z (first 
row), y (second row), and for the total intensity (third row). 
The excitation polarization of the external field is along z. For 
a scatterer with an isotropic Raman tensor, for example, the 
integrated intensity on the surface of the disk for \Ey^ divided 
by that for |-Ez|^ defines the depolarization ratio p. By visual 
inspection, it is evident that the depolarization ratio will be 
smaller at 312 and 364 nm than at 339 nm, where the disk has 
its main resonance. Colors mean different intensities for the 
different columns; we provide the maximum field intensity for 
each wavelength at the top of each column. Red (blue) means 
high (low) intensity in a linear color scale. See the text for 
further details. 
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FIG. 3: Depolarization ratios (top) and total SERS intensities 
(bottom) as a function of wavelength calculated for an ellipse 
for some of the models presented in Fig. As opposed to the 
case of a disk in Fig. Q an ellipse has two main shape-related 
plasmon resonances in the SERS intensity profile. The Tan — 
B model for the orientation of the molecule and the scattering 
mechanism gives depolarization ratios which go well above 1 
across the visible range. 
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FIG. 4: Comparison between the depolarization ratios (top) 
and the total SERS intensities (bottom) calculated for a 3D 
sphere with those of a 2D disk. The curves with symbols 
(lines) are the 3D (2D) cases for models Iso — B, Rand — B, 
and Rand — C, respectively. The depolarization ratios to- 
gether with the SERS intensities in this figure complete the 
qualitative/quantitative comparison of the effect of dimen- 
sionality on these quantities. 
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FIG. 5: Depolarization ratio (top) and SERS intensities (bot- 
tom) calculated using exact Mie theory for 3D spheres of var- 
ious radii a from 1 to 60 nm. Thicker (thinner) lines represent 
smaller (larger) radii. The results are shown for the I so — B 
case. 
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FIG. 6: Depolarization ratios (top) and SERS intensities (bot- 
tom) for Iso — B, Rand — B, and Rand — C, for a 2D dimer 
with separation d = 0.1a. Note that the depolarization ratio 
is not affected by all the structures appearing as resonances 
in the SERS enhancement. 
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FIG. 7: Depolarization ratios p (top) and SERS intensities 
(bottom) for the case Rand — C for different model geome- 
tries. This figure summarizes the various effects on p of shape, 
dimensionality, retardation, and coupled resonances. 



FIG. 8; Two scanning electron microscope images of Ag- 
colloids used in the experiments reported here and prepared 
by the method of Ref . |25| . The images are at different mag- 
nifications and give an idea of the typical spread in sizes (a) 
and the details of individual particles (b). 
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FIG. 9: Experimental depolarization ratios for 3 different 
laser excitations: 457, 514, and 633 nm, for a colloid sam- 
ple with (empty circles) and without (solid circles) KCl. The 
analyte is the benzotriazole dye BTZ2. These data show a 
dispersion in laser wavelength as well as a change depending 
on the interacting state of the colloids (KCl concentration). 



